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The phase transition of hadronic to quark matter and the boundaries of the mixed hadron-quark 
coexistence phase are studied within the two Equation of State (EoS) model. The relativistic 
effective mean field approach with constant and density dependent meson-nucleon couplings is used 
to describe hadronic matter, and the MIT Bag model is adopted to describe quark matter. The 
boundaries of the mixed phase for different Bag constants are obtained solving the Gibbs equations. 

We notice that the dependence on the Bag parameter of the critical temperatures (at zero chemical 
potential) can be well reproduced by a fermion ultrarelativistic quark gas model, without contribu- 
tion from the hadron part. At variance the critical chemical potentials (at zero temperature) are 
very sensitive to the EoS of the hadron sector. Hence the study of the hadronic EoS is much more 
relevant for the determination of the transition to the quark-gluon-plasma at finite baryon density 
and low-T. Moreover in the low temperature and finite chemical potential region no solutions of the 
Gibbs conditions are existing for small Bag constant values, B < (135 MeV)'^. 

Isospin effects in asymmetric matter appear relevant in the high chemical potential regions at 
lower temperatures, of interest for the inner core properties of neutron stars and for heavy ion 
collisions at intermediate energies. 
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1. Introduction 



It has been suggested that a phase transition would take place from hadronic matter to deconfined 
quark-gluon matter at sufSciently high density and/or high temperature. There has been considerable 
interest in relativistic heavy-ion collisions which could offer the possibility of producing a hot and 
dense matter and/or plasma of deconfined quarks and gluons. Up to now, experimental data on 
the phase transition have been extracted from ultrarelativistic collisions of almost isospin-symmetric 
nuclei, having a proton fraction Z/A ^ 0.4 - 0.5, in a region of very small chemical potentials. These 
experimental data from ultrarelativistic collisions could provide the possibility of studying the phase 
transition and testing effective QCD models. 

The interest in the transition at high baryon and isospin densities has been recently growing from 
the possibility of using new heavy ion facilities at intermediate energies and for the interest in the 
properties of the inner core of neutron stars where the transition to deconfined quark matter can likely 
occur. Now the main problem is that we do not have fully reliable effective QCD thoeries able to 
describe the two phases so the main approach has been based of a two-EoS model with the Gibbs 
conditions. In fact this scheme has been widely used to make predictions on the phase transition in 
the interior of neutron stars (e.g., se the recent Q and refs. therein). We remark that, even in 
the two-EoS approach, only a few papers have studied the phase diagram of hadron-quark transitions 
at high baryon density in connection to the phenomenology of heavy-ion collision in the ten A GeV 
range (intermediate energies) 

In Ref. 0] the phase transition from hadron to quark matter has been firstly analyzed also for isospin 
asymmetric matter. In this work, a Relativistic Mean Field (RMF) model, involving the interaction 
of baryons with isoscalar scalar and vector fields and with the isovector p meson and pion field, was 
used for hadronic matter and a MIT-Bag model involving massless u and d quarks was adopted for 
quark matter. More recently, Refs.0, [J, the RMF approach was extended to the isovector-scalar 
(5-field to study symmetry energy effects on the possible formation of a mixed hadron-quark phase at 
high baryon density during intermediate energy collisions between neutron-rich heavy ions. 

We remind that the nonlinear Walecka model (NLWM) 0] based on the RMF effective theory, 
has been extensively applied to study the properties of nuclear matter and neutron star, stable nuclei 
and then extended to the drip-line regions. In the last years some authors 10l4l4| have stressed the 



importance of including the (5(ao(980)) field in hadronic effective field theories for asymmetric nuclear 
matter. The role of the S meson in isospin channels appears relevant at high density regions 10l-|l4| 
and so of great interest in nuclear astrophysics. 

In order to describe the medium dependence of nuclear interactions, a density dependent relativistic 



hadron (DDRH) field theory has been also suggested |15l4l7| . The density dependent meson- nucleon 
couplings are based on microscopic Dirac-Brueckner (DB) calculations 16|, ll8| and adjusted to re- 



produce some nuclear matter and finite nuclei properties |15l - ll7l |. Recently the density dependent 
coupling models have been applied to the neutron stars 0. 

In this paper, we use the constant coupling (NL-RMF) scheme and the DDRH model for hadronic 
matter and the MIT bag model for quark matter. We study the phase transition of hadronic to quark 
matter and in particular the dependence of the boundaries of the mixed hadron-quark coexistence 
phase on the choice of the Bag pressure. Our study leads to the appearance of a Critical-End-Point 
of the mixed phase for low values of the Bag parameter. Moreover recently a similar effect due to a 



density variation of an effective Bag constant has been observed using a Nambu-Jona Lasinio (NJL) 
model, with dynamical varying masses for the quark phase [2o| . 

We also want to see the effects of the different hadronic models on boundaries of the mixed phase. 
We find that the critical temperature at low chemical potential is almost not affected by variation of 
the hadron EoS and the Bag-constant dependence can be well reproduced by a naive ultrarelativistic 
massless quark gas model. Conversely the critical chemical potential at vanishing temperature is 
very sensitive to the EoS of the hadronic part. Furthermore the mixed phase structure at high 
baryon density is affected by the hadron EoS, especially in isospin asymmetric matter. This has been 
observed when the isovector-scalar S meson is included and when some density dependence of the 
isovector meson couplings are introduced. 

The present paper is organized as follows. In order to give more emphasis on the new obtained 
results we report in the Appendices all the details about the effective Equation of State adopted in the 
two sectors: in Appendix A both models used for hadronic matter are described while in Appendix B 
we introduce the MIT-Bag model at finite temperature for quark matter. 

The Gibbs phase transition conditions are presented in Sect. 2. The results are shown and discussed 
in Sect. 3. Sect. 4 is devoted to a detailed analysis of the Bag constant dependence of critical end points 
in temperature and chemical potential. Finally some conclusions are drawn in Sect. 5. 



2. The phase transition 



According to the Gibbs conditions for phase transition, the temperatures, the chemical potentials, 
and pressures of hadronic matter have to be identical to that of quark matter inside the mixed phase. 
Moreover we must require the conservation of the total baryon and isospin densities : 

Aif(T,pf,pf) = /.g(r,pg,p?) , (f) 

M3^(T,pf,P3«) = ^3«(r,pg,p?), (2) 

P"{T,p^,p^)^pQ{T,p%p^) , (3) 

Ps = (l-X)pf + Xpg , (4) 

pJ-(l-x)pf + XP?, (5) 

where x is the fraction of quark matter in the mixed phase. 

The densities and chemical potentials for hadronic matter are defined as 

Pb = Pp + Pn , P3 = Pp~ Pn , (6) 

P-B ^ ^'^P-P + P-n) , P3^2^Pp-Pn), (7) 

and consistently for quark matter we have 

PB = ^iPu + Pd), p^=Pu-Pd, (8) 

3 I 

P%^2^Pu + Pd), P3 ^ ^{Pu - Pd) , (9) 



The asymmetry parameters for hadronic and quark matter are defined, respectively, 

H Pn — Pp Q ^Pd — Pu /, „N 

a = , — i , (10) 

Pn + Pp Pd+ Pu 

and the total asymmetry is 

«^ = -pI/pI- (11) 

Nucleon and quark chemical potentials, as well as the pressures in the two phases at the fixed 
asymmetry and a fixed temperature, can be obtained directly from respective EoS, as shown in detail 
in the Appendices A and B.. 



In the hadron sector we will use the Non-Linear Relativistic Mean Field models, 10|, [l^ |2l|, with 



different structure of the isovector part, already tested to describe the isospin de pen dence of collective 



flows and meson production for heavy ion collisions at intermediate energies, [l3j, |23| . We will refer 
to these different Iso-Lagrangians as: i) NL, where no isovector meson is included and the symmetry 
term is only given by the kinetic Fermi contribution, ii) NLp when the interaction contribution of an 
isovector- vector meson is considered and finally iii) NLpS where also the contribution of an isovector- 
scalar meson is accounted for. See details in Appendix Al. 

We are well aware that there are several uncertainties on the stiffness of the symmetry energy at 
high baryon density, mainly due to the lack of suitable data, see the reviews Therefore we 

will show also results with effective hadron interactions based on RMF models with density dependent 
meson-nucleon couplings {DDRH forces. Appendix A2) that present much softer symmetry terms at 
high baryon density. 

As already mentioned, in the quark phase we use the MIT-Bag Model, Appendix B, where the 
symmetry term is only given by the Fermi contribution. An important part of this work is the study 
of the Bag-constant dependence of the transition. 



3. Results: Properties of the Coexistence Zone 

The EOS for hadronic and quark matter is used to calculate the baryon and the quark chemical 
potentials, as well as the pressures in the two phases at the fixed isospin asymmetry and temperature. 
We study the phase coexistence region in the (T, pB, ps) space and even the properties inside the 
mixed phase [5|. For a fixed value of the total asymmetry and temperature T, we find solutions of 
the Gibbs conditions for both cases of x=0.0 and 1.0. In this way the boundary of the mixed phase 
region in the (T, pb) plane can be obtained. Finally we can get the full phase transition diagram by 
repeating the procedure. 

In our calculations we have used m„ = = 5.5 MeV. We will only consider the two-flavor case 
{q — u, d) in the bag model and the interactions between quarks are neglected as a first step. In fact 
in this work we are mainly interested in the study of the Bag Pressure dependence of the results. In 
order to obtain the binodal boundaries and the phase diagram, an iterative minimization procedure is 
adopted in the search for the solutions of the Gibbs conditions. In Figll]we present the phase diagram 
of hadronic-quark matter coexistence in (T, p, p^) space with 5^/^=160 MeV in the NLp and the 
NLpS models. In order to see the effects of different RMF models on the phase diagram, we also use 




the DDRH — RMF model to describe the hadronic matter with the same bag constant to calculate 
the mixed phase coexistence diagram, the results are given in FiglJJ 

From Figsll]and[2] we notice that the main differences between NL — RMF and DDRH — RMF 
models are the boundary densities {p%/pi3 and p% / pa) for the quark and hadronic phases at zero 
temperature in the isospin asymmetric case. Both models give almost the same critical temperature 
Tc=130 MeV, which then is not related to hadronic EoS differences as well as to the isospin asymmetry 
of the matter. 

In order to know the effects of different values of the bag constant on the hadronic-quark matter 
coexistence in (T, p, p^) space, we use B^/^—\%Q MeV and the two RMF models to calculate the 
phase transition diagrams; the results are shown in FigslJ] and HI respectively. Like in the previous 
B^/^=im MeV case, we can see that the behaviors of the binodal boundaries for symmetric matter 
in the quark and hadronic phases in (T, p) space given by the both NL — RMF and DDRH — RMF 
models are roughly the same. 

At variance, for asymmetric matter interesting isospin effects are appearing on the critical points 
at high baryon density and low temperatures, directly related to the high density behavior of the 
symmetry energy in the hadron phase [sl. [2o|. see also the discussion in Appendix A. In fact we get 
that for both Bag-constant choices the boundaries of in (T, p) space for the hadronic phase 

given by DDRH p and DDRHpS models are different from that given by NLp and NLp5 models. 

From Figs. [H [21 [3] and HI we remark that some squeezing-reopening effects of the binodal surface 
appear. The squeezing-reopening effect is more evident in the _B^/^= 160 MeV case and this is due 
to the variation of the pressure in the quark phase, for a fixed chemical potential /i^, as shown in the 
following. 
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FIG. 3: Same as in FigUwith B^'''^=190 MeV in the NL - RMF model. 



In order to simplify the discussion we consider the case of symmetric matter, Q!"^=0.0, and only the 
hadronic NL interaction. First we define the transition chemical potential = In FiglS] 

we plot the transition phase diagrams for symmetric matter in the (T, /i^) space with B^/^= 160 and 
190 MeV, in the NL model. The inset is the corresponding phase diagram of the transition in the 
(T, p) space for symmetric matter with 5^/^= 160 MeV. In order to discuss the squeezing-reopening 





FIG. 5: T ~ jiB phase diagram in the NLp model with £^''''=160 and 190 MeV. The inset is the T — pB phase 
diagram in the NLp model with B^''* = 160 MeV 



effect both (T, /is) and (T, p) phase diagrams should be considered. 

If we look at the (T, p) inset in FiglH we note that three main regions exist : 

(1) The "opening" region at low baryon densities and high temperatures. We can see from FiglS] 
that the transition chemical potential /i^ is smaller than the effective nucleon mass and the effective 
chemical potentials ^* even smaller in this region. This indicates that E* > /i* for the hadronic part 
and the nucleon fermion distribution will be small. So the p|| will present a slow increase with hb 
{fi{k) ~ fi{k)). We note that the hadron pressure is mostly of thermal nature and with a positive 
antifermion contribution in this region. The in the quark part is always much larger than the 
current quark masses (to„ = md—5.5 MeV) and so we can have a fast increase of the p^. 

(2) The "re-opening" region at high baryon densities and low temperatures. In this region, /is is 
larger than the nucleon effective mass. This means that E* < /i* for the hadronic part, the nucleon 
fermion distribution will not be small and so the will show a fast increase with /xb, also because the 
antifermion negative contribution is of course reduced. The interaction part of the hadron pressure 
shall increase. In correspondence we need a fast increase also of the p^ in order to keep the = 
balance. 

(3) The plateau at intermediate region of baryon densities and temperatures. In this region we get 
Pb — M* (nucleon effective mass), thus we can have the squeezing of the mixed phase in the T — pB 
plane and a kind of plateau in the T — pB plane. Our study clearly shows that the squeezing- reopening 
effect on the mixed phase is very sensitive to the value of the MIT bag constant. 
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FIG. 6: P — pB phase diagram of hadronic-quark matter for symmetric matter (ar^O.O) at temperature 
T = 35 MeV , obtained by the NLp model and various choices of the MIT-Bag constants. 

Finally in Figl6]we present the pressure-baryon density phase diagram of the hadron-quark transition 
for symmetric matter [ax—O-Q) always in the NL model and for various choices of the MIT-Bag 



constant. The temperature is fixed at 35 MeV, i.e. in the region of the mixed phase squeezing. 
The Maxwell constructions are indicated by the connections between the solid line of the hadron 
pressure and the dashed, dot-dashed and dotted lines of the quark pressure. It can be seen that the 
squeezing-reopening effect is larger for lower the bag constant choices. In fact for small values of the 
bag constant we cannot have the solution in the intermediate, flattening, region. This means that the 
Gibbs conditions cannot be satisfied and we get something like a critical End-Point of the mixed phase. 
From the figure we see that this is first happening for a Bag constant B^^^ ~ 155 MeV. If we further 
decrease the Bag pressure we reach the limit of no solutions in general at high chemical potential, 
since the baryon pressure cannot match the quark pressure, and the mixed phase just disappears. 
This important result at low temperatures and high densities will be further discussed in detail in the 
next Section, in particular in connection to the Fig llll 
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FIG. 7: Baryon and quark densities as a function of the bag constant at T=0. 




4. Dependence of the binodal surface on the Bag constant 

We will focus our analysis on results for isospin symmetric matter. In FiglT] we plot the density 
boundaries of the binodal region at T=0 as a function of the bag constant B^^'^ in the interval 140-240 
MeV. The dotted line is the quark matter limit density p%/ po, the solid line is the hadronic matter 
limit density p^l Po- We clearly see that both limits are increasing with larger Bag values, as expected 
from the fact the the energy per particle in the quark phase is increasing. It can be also seen that 
Pb/Po increases more quickly with the increasing Bag constant than the p^ / pa- So the selection of 
the Bag constant value will influence the binodal surface and in general the flnal results of the phase 
transition. 

This discussion about Bag constant effects on the transition is more clear in the {T,pb) plane. 
Already from FiglS] (for symmetric matter), where results are reported for B^/'^= 160 and 190 MeV 
values, it is evident that the transition region is deeply affected by the choice of the Bag constant. 



in particular at the end points at zero temperature or zero chemical potential. We will separately 
analyze the two regions. 
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FIG. 9: Critical chemical potential for symmetric matter (aT=0.0), at zero temperature, as a function of the 
bag constant B^^*. 



In Fig|S]we report the critical temperatures for symmetric matter, at zero chemical potential, for 
different values of the MIT Bag constant in the 160-220 MeV interval, using the NL model for the 



hadronic sector. The square points are the calculated values while the solid line is given by the simple 
linear fit : Tc = O.SlSSi/". 

Similarly in the Figl9] we show the critical chemical potentials, at zero temperature, for different 
values of the MIT Bag constant in the same 160-220 MeV interval, always for symmetric matter in 
the NL model for the hadronic sector. The square points are the calculated values while the solid line 
corresponds, even in this case, to a rather simple linear fit : /ic = a + with a = GSOMeF and 

h — 5.7. It is interesting to analyze such almost linear behavior of the two end-points as a function of 
the Bag values, see the following subsection. 

Finally, as already noted, we must remark that no /ic ^ solutions are existing at zero temperature 
for very small Bag constants, B < (135 MeV)'^. This is due to the fact that at low temperature and 
final chemical potential the hadron pressure, cannot compensate the quark pressure in the coexistence 
region, as we can also see from Fig[TT] Of course this effect is depending on the EoS in the hadron 
sector. This point will also be discussed in the next subsection. 



4.1 Ultrarelativistic gas model 

For the quark sector we can use as a reference the naive fermion ultrarelativistic gas model, con- 
sidering only free quarks in the Bag. Consistently with the results obtained before, no gluons are 



included. The pressure is simply given by 25l. |26| 



3 ^120 + 4^r^ ^ Stt^^t' ^ ^ ' 



where gq = gq = 12 is the spin-color-flavor degeneracy factor. The condition P"^ — 
will give the critical line of the transition to the quark deconfined matter. Neglecting the hadron 
contribution — condition), from the Eq. (fT2|) we can easily get the two critical end points, at 

zero chemical potential and at zero temperature, as a function of the Bag parameter. 

At /i = we get a critical temperature 

T, = [J^r-li/45i/4 ^ 0.81^1/4 MeV (13) 

in a very good agreement with the results obtained from the Gibbs conditions, solid line in FiglS] 
The hadron EoS is not contributing to the critical temperature at zero chemical potential, as we can 
also clearly see from the Figs[T] and [2j The reason is that at /i = only thermal vacuum excitations 
contribute to the pressure, ruled essentially by the particle degrees of freedom, which are much larger 



for the quark phase 26|. 



This is nicely confirmed by the Fig HO] where we plot the temperature evolution of the pressure in 
the two phases at /i = 0. We see that the hadron pressure is almost negligible up to the crossing 
points (the Tc obtained from the full Gibbs conditions). We can estimate an increase at much larger 
temperatures, due to the onset of anti-hadron contributions, also related to a decrease of the effective 
nuclcon masses.. 

Of course at fi — one can argue that the nucleon mass is too high to be thermally excited and in 
fact it is known that the hadronic matter is now dominated by pions. It is then interesting to compare 
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FIG. 10: Symmetric matter (aT=0.0) at zero chemical potential: Hadron pressure (solid, NL choice) and 
Quark pressure Eq.dll}, for various Bag constants (dashed, dot-dashed) as a function of temperature. 
The dotted line gives the pressure of a massless pion gas, see text. The crossing points correspond to the 
critical temperatures. 



the hadron pressure obtained from the RMF approach with the results of a pure massless pion gas 
pressure P^, dotted line of FigHUl given by |26,] 

^90 

where = 3 is the pion degeneracy. The critical temperatures are almost not affected. A more 
refined model would be to include all the mesonic resonances with their masses. However this does 



not significantly change the hadronic pressure up to T ~ 150 — 160 MeV [27|, |28|. We nicely see that 



the pure thermal hadron pressure is only ruled by the particle degrees of freedom, as already noted. 

At variance if we perform the same analysis at T = in order to get the B-dependence of the critical 
quark chemical potentials we get 



= [27r2]i/4si/4 ^ 2.11B1/4 MeV (14) 
which corresponds to baryon chemical potential, Ea. (|B7p . 



= 2 (/^c + Mc) = = MeV (15) 

which is in large disagreement with the results obtained from the Gibbs conditions, solid line of 
FigEl corresponding to a linear fit /x^ = a + bB^^^, with a = 680 MeV and b = 5.7, we note two 



differences: i) The B'^/^ slope is below the 6.33 value expected by the pure quark gas model; ii) We 
need a constant term a = 680 MeV. Both points are related to density dependence of the hadron 
interaction, that now cannot be neglected. The different slope comes from interaction contributions 
to the pressure in the hadron phase. The added constant is due to the presence of the nucleon rest 
mass in the hadron chemical potential. It is difficult to work out an analytical derivation of these 
effects since the relation between pressure and chemical potential in the hadron sector is not simple 
and naive Fermi gas models are largely overestimating the pressures in P — plots [29|. We can 
only state that at zero temperature and large chemical potentials the critical points are very sensitive 
to the hadron contribution. 
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FIG. 11: Symmetric matter (q:t=0.0) at zero temperature: Hadron pressure (solid, NL choice) and Quark 
pressure P'^ for various Bag constants (dashed, dot-dashed) as a function of baryon chemical potential. We 
show also P'^ = line corresponding to the case without hadron contribution. The crossing points (solid 
with the hadron part, open without the hadron contribution) correspond to the critical chemical potentials. 

This can be easily seen from FiglTll where we show the fully calculated pressures P*^ and vs. 
the baryon chemical potential at zero temperature. We show also the points on the line P*^ = 
corresponding to the deconfinement conditions without hadron contributions. We clearly notice the 
difference in the crossing points with and without the hadron part. The critical chemical potentials 
get rather larger values (of about a 50%) with the hadron contribution and also the B-dependence 
looks different. 

Finally we can work out a simple consistency check of the results shown in the Figs. [7] and IHl using 
the zero temperature relationship between quark chemical potentials and densities pq — /i^/tt^ . Since 
Pb = ^{Pd + Pu) |pg we finaly get 

Mg = 3M, = 3(^)V3(^)i/3(^e)MeF 



Using the critical quark chemical potential parametrization of FiglQl /ic = a + bB^^^, with a = 
680 MeV and b = 5.7, we nicely get the B^^^ dependence of the quark baryon density limit line of 
the coexistence region shown in FiglT] 



5. Concluding remarks 

We use the two kinds of RMF models for hadronic matter, with constant and density dependent 
couplings, and the MIT-Bag model for the quark matter to study the boundaries of the mixed coexis- 
tence phase in the transition from hadronic to quark matter. Binodal surfaces and critical end-points 
for different isospin asymmetries are obtained using different hadronic models and various MIT Bag 
constants. The Bag Pressure dependence of the critical points is very interesting since we clearly show 
the difference between two regions: 

• High Temperature and Low Baryon Chemical Potential. 

The hadron EoS is not important: we get the same Tc for any hadron EoS, even including a 
pure massless pion gas contribution. Isospin effects in asymmetric matter are also negligible. In 
conclusion the critical temperature at zero chemical potential, when existing, only depends on 
the value of the Bag constant. In fact we show that if small Bag values are used the binodal 
surface is shrinking at low densities and finite temperatures and eventually no continous solutions 
down to ^ = are found for B < (150 MeV)^. 

We note that a Tc comparison with recent lattice QCD calculations with physical quark masses 



28| would indicate an effective Bag constant, at low fi, around B ~ (180 — 190 MeV) 
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High Baryon Chemical Potential and Low Temperature. 

Now we have just the opposite evidence: the Critical Chemical Potential is clearly dependent 
on the interaction in the hadron phase. Actually we even show that if small Bag values are 
used the binodal surface is shrinking at high densities and finite temperatures and eventually 
no /i 7^ solutions are found for B < (135 MeV)'^. 

Moreover for asymmetric matter all that translates into relevant isospin effects, very sensitive 
to different symmetry energy terms in the Hadron EoS. These results are of interest for all 
nuclear systems at high baryon and isospin density and moderate temperature, like neutron star 
formation and inner core structure and/or heavy ion collisions at intermediate energies. 
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Appendix A: Equation of state for hadronic matter at finite temperature 



Al. Nonlinear relativistic mean field model with constant couplings 

The Lagrangian density, with the isovector scalar S field, used in this work is 

+ l{d^6- d^S - m^/^) - \f^.F'^-' - \g^M^\ (Al) 

where </> is the (/i-meson field, is the w-meson field, 6^ is p meson field, 5 is the isovector scalar field 
of the ^-meson. = — d^^^^, G^^v = Q^Jiv — djb^^ and the J7(0) is a nonlinear potential of a 

meson : U{(j)) ^ ifl^^ + jbcf)'^- 

The EOS for nuclear matter with the isovector scalar field at finite temperature in the mean-field 
approximation (MFA) is given by 

^37 



i—n,p 



2E / + /<(*)) + 5"'i*' + fW 



-n,p 



where E* = \Jk'^ + M*'^. The nucleon effective masses are defined as 

Mi* = M — g^cj) =F gsSs (— proton, + neutron). (A4) 
The field equations in the relativistic mean field (RMF) approach are 

(t> = —^4>'' - ^<t>^ + ^{Psp + Pan) , (A5) 

mi mi mi 



wo = , (A6) 
mi 



h = —{Psp - Pan) , (A8) 

m^ 

with the baryon density p = = pp + p^, p3 = p" = pp - pn and ps3 = Psp - Psn, Psp and Psn are 
the scalar densities for proton and neutron, respectively. The fi{k) and ,fi{k) in Eqs.(A2)-(A3) are 
the fermion and antifermion distribution functions for protons and neutrons (i = p,n): 



and 



^'(^^-i + exp{[£;>(fc) + Mt]/r}- (^''^^ 

where the effective chemical potential ^i* is determined by the nucleon density pi 



(27r) 



:mk)-m)), (All) 



and the fi* is related to the chemical potential /ii in terms of the vector meson mean fields by the 
equation 

p* — fii ~ gujLOQ =F gpbo (~ proton, + neutron), (A12) 

where p,i are the thcrmodynamical chemical potentials fii = de/dpi. The chemical potentials for 
proton and neutron are given by, respectively 

tJ-P^ y-p^ jP H tP3 

Mn = + 4-^ - ^^3. (A13) 

The proton and neutron chemical potentials can be denoted in terms of the baryon and isospin 
chemical potentials by the equations 

/^p = AiB + M3, Mn = Ms-M3- (A14) 

The scalar density ps is given by 

r Aih ]\,T* 



i—n,p 



A 1.1 Parameters 

The isovector coupling constants, p-field and p + S cases, are fixed from the symmetry energy 
at saturation and from Dirac-Brueckner estimations, see the detailed discussions in refs. IC 



The coupling constants, fi = gf/mf, i = a,u!,p,6, and the two parameters of the a self-interacting 
terms : A = a/g^ and B = b/g"^ are reported in Table 1. The corresponding properties of nuclear 
matter are listed in Table 2. Here the energy per nucleon is defined E/A = e/p — M. 



Table 1. Parameter set. 



Parameter Set 


NLp 


NLp(5 




10.329 


10.329 




5.423 


5.423 




0.95 


3.150 




0.00 


2.500 




0.033 


0.033 


B 


-0.0048 


-0.0048 



Table 2. Saturation properties of nuclear matter. 



Po (/m ^) 


0.16 


E/A {MeV) 


-16.0 


K {MeV) 


240.0 


E,ym (MeV) 


31.3 


M*/M 


0.75 



A2 Density Dependent coupling model 



The Lagrangian density, with 6 meson, now reads 



+ ^-{8^5 ■ d'^S- m^J^) _ _ Ig^.G'^^ (A16) 

where (p is the (^-meson field, the w-meson field, the p meson field and 6 the isovector scalar 
field of the i5-meson, respectively. F^^, = df^ujiy — dj^uj^ and G^^ = d^jb^, — 9^6^. 

The most important difference to conventional RMF theory is the contribution from the rearrange- 
ment self-cncrgics to the DDRH baryon field equation. The meson- nucleon couplings g^^ g^j, gp and 
gs are assumed to be vertex functions of Lorentz-scalar bilinear forms of the nucleon field operators. 
In most applications of DDRH theory, these couplings are chosen as functions of the vector density 

The variational derivative of Lagrangian density Eq. (A16) is 



dc_dc^dcdi 

d^p d'lp dp d-ip 

The density dependence of the vertex functions go-, goj, gp and gs produces the rearrangement 
contribution 



(A18) 



and 



e;,^ - u^, (A19) 
where u^=^ is a four-velocity with m^u^=1. The Dhac equation takes the conventional form 

[7^(za^-E^)-(Af-S^)]^ = 0, (A20) 

where = Sw'^/j + Sp^^' ^p ~ and = 9a4> + 95^- S. The field equatfon of baryon in the mean-field 
approximation (MFA) is 



(A21) 



with 



1 9cr 5ct / I \ 

= —Ps = -^[Psp + Ps,i), 

mi mi mi 

bo = ^ < -ipl^Ts^jj >= -^p3 = -^(pp - p„), 

r 9S ^ J , ^ 95 _ gs , . 
mi mi mi 



ig nig nig 

vfl. fdga.ga 2./9gs gs 2 ^dg^ g^, 2 (dgp.gp 2 /Aoo^ 

where pi(i = p,n) and psi are baryon and scalar nucleon densities, respectively. 
The EoS at finite temperature in the MFA is given by 

i=n,p V y a- uj p d 

>^ 2 /■_d3fc fc2 1 5^ 2_L 1 2_L 1 -^P 2 1 55 2 v« tAOA\ 



where Ei* = y fc^ + A/^*^ i=p,n. 

The proton and neutron chemical potentials in DDRH model are 



/^p — Pp 


, gl 

\ 2"/'" 

mi 


mf, 


^0 


Pn = Pn ^ 


9l 

' —P ' 
mi 


9l 

- —P3 - 


^0 ■ 



(A25) 

and the baryon and isospin chemical potentials in the hadron sector consistently defined like before, 
Eq.dJTll. 



A2.1 Parameters of the meson-nucleon coupling density dependence 



The parameters of the model include nucleon mass M — 939MeV, the masses of the mesons mg., 
m,.; , nip, ms and the density dependent meson-nucleon couplings. The parametrization was proposed 
as : 



t{p) = 9i{Po)h{x), for i = a,uj,p,d, 



(A26) 



with 



, l + bi{x + di)^ 



(A27) 



where x — pj po and po is the saturation density. For p and d mesons, the following parametrization 
was proposed [3l| : 



fi{x) = a.iexp[~bi{x - 1)] - Ci(x - di), for p,S. (A28) 

The parametrization form and parameters are taken from ref. 17 1 for a, uj mesons and ref. [si^ for p, 
S mesons, respectively. All parameters are listed in Table 3. 

Table 3. Parameters of the DDRH model. 



Model 


TW [17] 


DDRHp [28] 


DDRHp^ [28] 


Meson 


a 




P 


P 


5 


rrii (MeV) 


550 


783 


770 


770 


980 




10.72854 


13.29015 


3.587 


6.476 


7.58963 




1.365469 


1.402488 


0.095268 


0.095268 


0.01984 


h 


0.226061 


0.172577 


2.171 


2.171 


3.4732 


Ci 


0.409704 


0.344293 


0.05336 


0.05336 


-0.0908 


di 


0.901995 


0.983955 


17.8431 


17.8431 


-9.811 



The corresponding density dependence of the couplings is presented in the FigfT2l 



A2.2 Hadronic EoS without and with density dependent couplings 

In order to better understand the effects of the density dependence of the nucleon-meson couplings 
on the hadron-quark phase transition we look here in more detail at the differences on the hadron 
EoS, NL vs DDRH approaches. In Figs. [T31 [HI and [T51 we show in the order the results for the 
energy per particle, the pressure and the baryon chemical potential in the case of isospin symmetric 
case at zero temperature. With density dependent couplings we see a more repulsive hadronic matter 
and so in general we would expect an "earlier" transition to the quark phase with increasing density. 
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FIG. 12: Density dependence of the meson-nucleon couplings in the used DDRH interactions. The bottom gp 
dashed hne corresponds to the DDRHp case, without the S meson. 
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FIG. 13: Energy per nucleon as a function of the baryon density for symmetric matter at T=0. Solid line: 
NL results. Dashed line: DDRH. 

In any case for symmetric matter the effect will not be too large. In fact we note that the DDRH 
pressure is not much larger that the NL values. The same is happening for the baryon chemical 
potentials. This is mostly due to the increase with density of the rearrangement corrections S^/O, 
Eq. (|A24|) , as shown in the inset of FigfTil 

At variance the effect of density dependent couplings will be noticeably larger in the case of isospin 
asymmetric matter. In fact at high baryon density the symmetry term will be much reduced due to 
the decrease of the p — nucleon coupling, see Fig[T2l This can be clearly observed from the FigfTHl 
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FIG. 14: Pressure as a function of the baryon density for symmetric matter at T=0. Solid line: NL results. 
Dashed line: DDRH. In the inset we show the density dependence of the rearrangement term correction to 
the DDRH pressure 
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FIG. 15: Baryon chemical potential as a function of the baryon density for symmetric matter at T=0. Solid 
line: NL results. Dashed line: DDRH. 




Appendix B: Quark matter equation of state 



The energy density and the pressure for the quark system are, respectively, given by the MIT Bag 
model 



32| 



e = 3 X 2 

q—u,d,s 



d^k 

(2^ 



(Bl) 



P 



3x2 

q—u,d.s 



(27r)3 



--{U + h)-B 



(B2) 



In 



where B denotes the Bag constant (the Bag pressure) chosen as a rather standard parameter 
fact a detailed discussion about the effects of different Bag constants on the hadron-quark transition 
will be the main topic of this paper, m^, q = u,d, are the free quark masses (5.5 MeV choice), and 
fq, fq represent the Fermi distribution functions for quarks and anti-quarks, respectively : 



1 



(B3) 



and 



l + eMiEq + fiq)/Ty 



(B4) 



where Eq = ^ + and /i^ are the chemical potentials for quarks and anti-quarks of type q. 
The quark number density is given by 



n,; =< qlqi >= 3x2 



/ 



d^fc 



i =: u,d, 



(B5) 



(27r)3 



where the chemical potential fiq can be obtained by Eq. (B5) for a given quark number density. The 
quark number density is Uq — Uu + = Pq and relation with the baryon density is 



The chemical potentials of quarks are related to the baryon and isospin chemical potentials 
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